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Abstract 

This work is an application of the second order gauge theory for the Lorentz group, where a 
description of the gravitational interaction is obtained which includes derivatives of the curvature. 
We analyze the form of the second field strenght, G = dF + fAF, in terms of geometrical vari- 
ables. All possible independent Lagrangians constructed with quadratic contractions of F and 
quadratic contractions of G are analyzed. The equations of motion for a particular Lagrangian, 
which is analogous to Podolsky's term of his Generalized Electrodynamics, are calculated. The 
static isotropic solution in the linear approximation was found, exhibiting the regular Newtonian 
behaviour at short distances as well as a meso-large distance modification. 
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I. INTRODUCTION 



Nowadays there are many proposals to modify gravitation in order to solve several prob- 
lems as the present day accelerated expansion of the universe , or to accommodate 
corrections of quantum nature which arise from the classical effective backreaction of quan- 
tum matter in a curved background j^. Effective action is widely used in quantum field 
theory as a powerful method of calculation. The Podolsky generalized electrodynamics, for 
instance, can be viewed as an effective description of quantum correction to the classical 
Maxwell Lagrangian 

For gravitation, usually higher orders terms are introduced by means of Lagrangian con- 
tributions quadratic in the Riemann tensor and their contractions j^. This is inspired by 
1-loop corrections in the Einstein-Hilbert action in the quantized weak field approximation, 
or in the equivalent Feynman construction of a spin-2 field on the flat Minkowski back- 
ground [g]. Besides this, at the quantum level, the S matrix for the Einstein theory is finite 
at one-loop level, but diverges at the two-loop order which motivates the introduction 
of derivative terms in the Riemann tensor for the action 8|. 

On the other hand, recently was proposed a second order construction of gauge theories 
based on Utiyama's approach 9|], which gives exactly the same correction terms as in the 
Podolsky electrodynamics, but now arising from the principle of local gauge invariance [l^ . 
Therefore, a connection between quantum corrections and gauge higher order terms in the 
action was conjectured, which was proved be fulfilled also for the effective Alekseev-Arbuzov- 
Baikov Lagrangian of the infrared regime of QCD 

Here, we analyse the gauge formulation of the gravitational field based on the framework 
of the second order gauge theory. The simplest gauge group is given by the Lorentz homoge- 
neous group in the context of a Riemannian description of the gravitational field. Since the 
gauge field is given in such case by the local spin connection, higher order in the gauge field 
involves naturally the derivative of the curvature tensor. In this sense, the actual higher 
order gravitational lagrangian should be constructed from invariants using the covariant 
derivative of the Riemann tensor instead of the usual quadratic terms in curvature. 

The relationship between the algebraic gauge description and the geometrical one is set- 
tled by means of the introduction of the tetrad field, and the construction of the covariant 
derivatives associated with the both symmetries: the local Lorentz and the global diffeo- 
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morphic coordinate transformations. We use Latin indexes, a,b,..., for the internal Lorentz 
group and greek indexes for the tangent space of the space-time manifold. 

The paper is structured as follows. In section [III we review some results relating gauge 
invariance and gravitation. The field strengths F and G of the second order treatment are 
introduced in section llllt where they are also written in their geometrical counterparts: the 
Riemann curvature tensor and its covariant derivative. 

Section IIVI deals with the possible quadratic invariants of the type and G^. All 
the possible contractions are studied and only the independent invariants are kept. This 
counting is made in the same spirit as the systematic selection of the independent Riemann 



monomials done in 121 and 



ISj. In the following section, section |Vl these invariants are 



shown to satisfy the identity which restricts the theories that may be called of the gauge 
type. 

Among all invariants, we select Lp = S^Rp^S^W^^, the Podolsky-like Lagrangian, 
for calculating the equation of motion of the gravitational field. This higher order gravity 
application is done in section IVIi For this Lagrangian, we calculate the static isotropic 
solution in the linear regime at section IVI Bl finding the regular Newtonian potential at 
short scales, but with a modified potential at intermediary scales. 

Final remarks are given in section IVIII 



II. GAUGE INTERACTION AND COVARIANCE 

In 1956 Utiyama ^ has shown how to implement a gauge description for gravitational 
interaction with matter fields (x) transforming according to 

sq^{x) = ^s^\x){j:,,)%q\ (1) 

as an implementation of the local invariance exigency of the action under continuous proper 
Lorentz transformations, which are characterized by the generators T^ab satisfying the oper- 
ation of a typical Lie group, 

1 

[Sab, Scd] = -fab, cd^ef , (2) 

where 
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are the structure constants obeying the Jacobi identity. = —e^"- are the parameters of 
the local transformation. The capital latin indexes are for the components of the matter 
field. 

It was clearly shown the need to introduce the compensating field cu^^ {x) transforming 
as a connection, 

= le''' (x) a{ ^.u^'^l + d,e^^ (x) . (3) 

To ensure the covariance under coordinate transformations it was necessary to define an 
space-time connection whose behaviour under infinitesimal diffeomorphisms is 

d5x- , d5x' d5x\^ dHx^ 

OV „^ = -—^L „^ 1 ^„. 1 , 



dx^ dx^" dxi'dx" 

The invariance of the theory implies that the compensating field must appear through 
the gauge covariant derivative 

D^Q^ = d,Q^ - lou^^ (E„,)l (4) 



I.e., 



5L'^g^ = ^£«^(E„6)lL'^g^, (5) 



and the space-time connection must appear through the space-time covariant derivative, 

s.Q'" ^ d,Q'^ + r^^g^^ + r;,g^" , (6) 

and the total covariant derivative: 

v.Q'" = d^Q"" - u;%,Q"' + r;„g- . (?) 

This total derivative must commute with the mapping to the tangent space of the manifold,^ 

Q'^ = h^Q'\ Q'^ = hiQ'" , (8) 
V^Q'" = h^V^Q'^ , (9) 



^ Note that the action of the total derivative on a tangent space field is defined by 
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where we have introduced the tetrad field h: 

Khj = Khj = S'j , g^i. = y]ijh\y^ , 7]ij = hi'^h/Qf,^ , h = -^/det^ = • 
The definition ([9]) imphes the absolute parallelism of the tetrad: 

V^/^i, = , (10) 
which can be solved for the compensating field, 

or for the space-time connection, 

= [D,hQ . (11) 

We will restrict our analysis to a symmetric space-time connection in order to approach 
the Riemannian description. The extension to the Riemann-Cartan case is quite natural, but 
would imply different types of invariants as admissible Lagrangians (see discussion bellow). 



III. GAUGE FIELD LAGRANGIAN 



The basic hypothesis we will assume is: the Lagrangian for the free gauge potential 
depends on the field, its first and second order derivatives and Lq = Lq {u!^j[, di^u'^j^, dpd^uj^ 



obeys local invariance under ([3]). This enable us to use the results presented elsewhere [10 1 
to construct a gauge formulation for higher order gravitation theories. 



A. The field strengths 



According to the work [10|, we can reexpress 

splitting it into a set of four hierarchical equations after substituting ([3]) and claiming the 
independence of the parameters e""^ and their derivatives. Three of these functional equations 
are used to conclude that 

Lo = Lo{F,G); ^^0, (12) 



where 



F^l = d.ujf - d^ujf - Vc^uj-^y^ + Vc,u;-\u;'\ (13) 



and 

G"',,. = DpF'^l = - mu^^F'^l^ + r//,^"/F^J^ . (14) 

The remaining hierarchical equation put in terms of the gauge fields F and G, 



Qpad-lbcgh^ pu ^ QQad J be gh'^ jipa — ^ ^ K^"^ ) 

imposes restrictions upon the functional form eventually chosen for Lq. Substituting the 
structure constants, this condition can be explicitly written as 

BT BT 

^ VlcgSl - r^.,^:] Ff^ + ^ [r^^,5l - r^,,5«] = . (16) 

B. Geometrical variables 

In this section we will show how to interpret all objects and condition of the previous 
sections in terms of a geometrical point of view. From ( fTOl) we read 

and therefore the field strength F is written as 

where we recognize the expression of the Riemann tensor [isj l. 

i.e., 

F^l = v'-K'-h^R^^^ ^ . (17) 

The easiest way to find the geometrical counterpart of G is to apply the geometrizing 
relations (I6]l8|) : 

ha'K-'G^^^ = K'h,-^D,F^l = 5,F^;^ ■ 

^ pa ~ "'a "'b ^ pa 
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(18) 



and use (fTTI) . We arrive at 

which is the most natural equation one would expect in view of the relation f|T^ between F 
and G. 

By means of the geometrical descriptions ( ITTll and ( ITSjl . we are able to find 



dLo _ dLo dRxp^^' _ dL, 







d i5pR^pa^ 



With these derivatives, the condition ( JTSj) for the gauge Lagrangian is put in the form 
dLo 



OR 



(19) 



+ 



dLn 



dLo 



5pR. 



'a pa. 



= 



This is a fundamental restriction upon the Lagragians tentatively proposed for the theory, 
and it is quite useful in order to choose a specific suitable invariant. 



IV. QUADRATIC LAGRANGIAN COUNTING 

Our goal here is to determine all possible independent quadratic Lagrangians constructed 
with the field strength tensors F and G considering their various symmetries. By quadratic 
Lagrangians we mean invariants of the type FF or GG, but not mixed terms like FG 
(obviously with the proper contraction of indices). We will also compute the linear case of 
the Einstein-Hilbert Lagrangian. 



A. First Order Invariants 

The symmetries to be considered in the construction of the invariants of the type FF 
are those inherited from F: skew-symmetry in each pair of indices: -F"^^, = —F^^u and 
F^-jlj^ = —F°'l^. Besides these, there is another which is unveiled by the geometrical form of 
F, eq. f lT7|) . namely 
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the familiar first Bianchi identity met in the context of the general relativity. 

Once algebra and space-time indices can be transformed into each other by means of a 
tetrad, we will consider a compact representation for F: 

Since the Lagrangians are all of the form with all allowed orders of contractions, it is 
always possible to rename dummy indices in such a way that the first F will keep its indices 
in alphabetic order. In the table below it follows all available permutations for the second 
F: 





Fix. a 


Fix. b 


Fix. c 


Fix. d 




{abed) 


{bacd) 


(cabd) 


(dabc) 


cyclic 


{acdh) 


ipcda) 


(cbda) 


{dbca) 




{adhc) 


{bdac) 


{edab) 


{deab) 




(abdc) 


{bade) 


{eadb) 


{daeb) 


non-cycl. 


(acbd) 


{bead) 


{ebad) 


{dbae) 




(adcb) 


{bdca) 


{cdba) 


(dcba) 



By means of a change in one pair of indices, one can see that the non-cyclic permuta- 
tions are all proportional to the cyclic ones. Considering only the cyclic permutations and 
changing two pairs of indices, the table is reduced to: 





Fix. a 


Fix. b 


Fix. c 


Fix. d 




{abed) 








cyclical 


{aedb) 


{beda) 








(adhc) 


(bdac) 


(cdab) 





The skew-symmetries of the first F (which has been taken in alphabetic order) leads one 
to restrict once more the possible contractions to the three quadratic invariants 

l[ = {abed) {abed) 

7f = {abed) {aedb) (21) 
l[ = {abed) {cdab) . 

We now analyze invariants constructed with one trace of F. The only non-null type of 
trace are those obtained by contracting one index of the first pair with one index of the 
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second pair, in view of the skew-symmetry of this object. All possibilities are proportional 

to 

TrF hJ^F^/^^hf^^ = i-ab-) or (oabo) . 
The quadratic invariants are given by, 

ifrF ^ (.^5.) (o^^o) (22) 

IJ^P = (.ah-) (obao) . 
Still, one can construct a linear invariant taking a double trace of F: 

jTrTrF ^ p^^^^hj^ ^ (■ O c) . 



B. Second Order Invariants 

Let us introduce a similar notation to the one used in the case of F, i.e., 

h/h/h/C^^^^ ^ [abode] , 
where we identify the following symmetries: 

(i) antisymmetry by permutation of indices in the first pair and last two of them, 

[abcde] = — [bacde] = — [abced] ; 

(a) Bianchi identity among the last three indices, 

[abcde] + [abdec] + [abecd] = . 

1. Invariants of GG kind 

The quadratic combinations are now in a larger amount than in the F"^ case. In fact, we 
have five tables like fl20l) . one to each letter labeling, since we can associate 

[abcde] = c (abde) . 
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Using the symmetries cited above, one finds that the 5! invariants are reduce to just 
two kinds: 



if = [abode] [abode] 



-'2 



[abode] [debao] . 

The detailed and cumbersome calculations are made in the appendix |Al 



2. Invariants Involving Traces 



There are three independent types of traces for G: 



^(1) 

abc 
(2) 



■a ■ bo] 



abc 
(3) 



hfG%^^hfhf = [-ab ■ o] 



Till = a'^G^p^h: ^ [ab ■ -c] . 
Again using symmetries (see appendix [X]) we arrive at: 



TrQ^ = [ab ■ o][-a ■ bo] 


TrQii = [ab ■ c] [bo ■ a] 


TrQ^ = [-ab ■ o][-o ■ ab] 


TrQi4^ = [ab ■ o] [ab ■ o] 


TtQq = [ab ■ -o] [-a ■ bo] 


TrQij = [ab ■ -o] [ab ■ -o] 


TrQio = [-ab ■ o] [-ba ■ o] 


TrQis = [ab ■ -o] [ao ■ -b] 



while for double traces we have: 



TrTrQi = [ o b ■ o][ o b ■ o] 

TrTrQ2 = ■ '°] ■ "o] 
TrTrQs = [ o b ■ o][ob ■ -o] . 



(23) 



(24) 



(25) 



C. Bianchi Identities 

As we already said, until now we have not used the first Bianchi identity: 

^<xxP/3 ^ j^xP-P + jiP-xfS ^ . (26) 

In geometrical variables, the cyclic property of G is translated to the second Bianchi 
identity: 
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These identities reduce the number of independent invariants, since F (x R and G oc SR. 



1. Reducing invariants 

Let us begin by invariants of form F^. The first three are (|2T|) : 

jF _ T3 TDcrXPK. 

As a consequence of the first Bianchi identity (l26l) and the skewsymmetries, the curvature 
tensor obey: 

^apXK Rxi^o-p • (^''') 

Then, 

jF _ TD T3Xi^<^P — E> U<^PXi^ — jF 

while for I2 one finds, 

jF _ T3 po-xpft _ _ (T3 _|_ p po-XPt — _ R PXO-PK _|_ p pxo-pt — —jF I jF 

which let us with only one invariant of this kind, . 

Now, we translate the trace-like invariants fl2^ in a geometrical form: 

= Rppu'Rr 

jTrF _ p p p upa 

-'2 ~ ^ppy -^o- 



Since the Ricci tensor R^^ = -^pmi/^ symmetric, we have in fact only one invariant, // 
RpvR^^ ■ 

At last, the only invariant of double traced form in F is: 

jTrTrF _ 



TrF 



^ Which is a consequence of the first Bianchi identity. 
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Analogously, in view of the Bianchi identities, only four invariants of the type remains 
(see appendix |A|) : 





TrTrg^ = SfRp^S^Rf'^ 







V. GAUGE INVARIANCE CONDITION 



With the invariants constructed above we collect seven types of Lagrangians for the 
gravitational field: 



Lagr. 


Inv. 


Gauge Form 


Geom. Form 


t{Ri) 
^0 


(^jTrTrFy^ 




i?" , n = 1,2 


t{R2) 


tTtF 
-'1 


bfia c 


RfiuR'^'^ 






rpab Tp 
fj.u-'- ab 


R n R'^l^f"' 


r(Gi) 


TrTrQ2 


^ Pa^cbfi 
^abp ^ fic 


6PRp^6^R>^^ 


r(G2) 




^ fiaa^cb 


5pR,^5^R''^ 


r(G3) 


TrQio 


riabe da 
^ aa^deb 


dfsR^J^R^f^ 


7-(G4) 

-^0 


I? 


r^ab r~i /''^-^ 
^ fiuX^ab 


^pR^p^J^R'^P^^ 



(28) 



We are considering Lagrangians only up to quadratic order in F and or G, which also 
include the linear invariant J'^^'^^^ = R and their square R^. Actually, one can observe that 
if any invariant fulfills the gauge invariance condition, then any of its power will also do it. 



since this condition is linear in the derivatives and For instance 



dLo 



dF 



r 



OF 



dG 

dF 



Therefore, 



9/ 



[...]F = 



dLo 
dF 



[...]F = 



and the same follows for G. 

Using the skewsymmetry z/ ^ 7 of the equation f|T9|) and the symmetry properties of the 
Riemann tensor, one can easily verify that all Lagrangians densities listed in ( 128|) accomplish 
the gauge invariance condition. Then, any function of these invariants expressible in a Taylor 
series also will fulfill the gauge invariance condition. 
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VI. EQUATIONS OF MOTION 



Here we will concentrate our attention on the effect of the term 

on a gravitational theory based on the Einstein-Hilbert action plus the L^^^^ term. This 
Lagrangian density is equivalent, by the Bianchi identity, to the form ^h6''Rp^6fj,R^^, which 
is clearly analogous to Podolsky's second order term for Electrodynamics {Lpodoisky oc 
d^Fp^dpF^'^). The choice of the particular Lagrangian Lq*^^^ is mainly motivated by this 
analogy. Besides this, the i^o*^^^ term also can be viewed as a kind of kinetic term for the 
scalar curvature, what approximate such description to the usual scalar fields. Moreover, 
this scalar is, up a surface term, present in the Schwinger-DeWitt renormalized effective ac- 

n 

tion for an scalar field on a curved background pj[. Therefore, the field theory constructed 
over basis on Lq can be considered an effective gravitational theory. 
Taking a functional variation of the tetrad field, one finds: 

h = ^ , 6h= hg'^'Sgxu = hg^^ h\r,ab5h\ 

and 



54^1) = -dp {hdPRSR) - -5Rdp {hdPR) + 



-g^'dPRdpR-g'^'^gP^dpRdpR 



On calculating the equations of motion, we must give special attention to the last term 
involving 

5R = -2R,pg^^^g^^h\7^,,5h\+ 



h 



h[g^^5Tl,-g''-5T^^^ 



which will include several integration by parts. After these integrations and some cumber- 
some calculations, one finds: 



SxS, [OR] + ^6xR6,R - Rx.OR - 9xu0 [OR] - ^gxJ'RSpR 
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where 



+ -h {g^-'g-^ - g'^'^g^^) 6^ [OR] {S'j\6\ + 6%5\5\ - 6%5\5\) 6g,, , 



and 



= Sp6^ 

is the Laplace-Beltrami operator on the Riemannian space. 

Therefore, the second order contribution to the equation of motion will be 

Hi = h'^6x6, [OR] + h'^6xRS,R - R\OR - h\5p5^ [OR] - hi\5PR5pR . (29) 

Furthermore, if we include the usual first order Einstein-Hilbert and a matter Lagrangian 
densities, 



St — I d"'X ( — — Lq + hCmatter] 

J \ X ) 
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the field equations become 

G\^i5E\ = xT\ , (30) 



or in a geometrical form, 

[OR) + ^SxRS.R - Rx,OR - gxuO [OR) - ^gxJ'RSpR = xTxu 
where G'',, is the Einstein tensor and 



R\u — -^gxuR + P 



Ti 2 (5 {hCYnatter ) 



h dg^" 

is the energy-momentum tensor of the matter fields written in terms of the metric field. 



ll| , one could expect that the higher order 



By analogy to the Alekseev-Arbuzov-Baikov 
terms, which can be until sixth derivative order, would be related to infrared corrections to 
General Relativity, giving sensible physical effects at large scales. 

A. Covariant Conservation of Txv 

Taking the covariant divergence of (l30l) . we have 
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Now, from the first order case, we know that 

Applying the divergence to the equation (!29|) . one finds 

d'H,^ = 5^5MR - g.J'^O [OR] + ^6,R5''SaR+ 
+^5''5,R5aR - S'^R.aOR - RuJ'^OR - \g.J" {5'RSpR) . 
Using the commutation relation 

[6,, 6a] A^ = RJ^A^ 
and the second Bianchi identity, we arrive at 

S^H^a = Ra^S^OR - RuJ^OR + ^S^RS'^Sc.R - h^RSJ^R = 
Then, the covariant conservation of T^jy is established: 

{Gp, + f3Hp,) = =^ 6^T^, = , 
as expected from the coordinate invariance of the Lagrangian density. 

B. Static Isotropic Solution 

In the case of an static isotropic metric, 

els'" = e-'^'-^dt^ - e^^'Ur^ - r^dO^ - sin^ Odd)'' 



in the vacuum, the equations of motion fl30|) are reduced, in the linear approximation, to 
the following coupled linear equations: 

+ 

\ 'j'^ ^ 'y lyO ijnO /y\) j 

I \ r r-^ r J 
( ^(6) + 3^(5) _ 1^^(4) + - -A(^) + 4a(^) + 8^ - 48^ + 96^ - 96 = o . 

\ ry ry A /yO /y ry A ryO ry^ ryi) ryK) j 
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To solve this system, we use the Probenius method, based on a series expansion: 



s+n 



Prom the first terms in the series, we find s = — 1, and the recursion relations above: 

i/„ (n - 2) 



4(3 {n + 4) (n + 2) {n + 1) (n - i) ' 



-'n+4 



2/3 (n + 4) (n + 3) (n + 2) (n - i) 
such that the solution can be written as 

3 / oo 



m=0 



n=0 



^ 1 ^ (4n + m - 2) (4n + 3 + m) 



2 (4n + 1 + m) 



m=0 m=0 n=0 

where z/^ and A^, with m e {0,1,2,3}, are the integration constants specified by the 
boundary conditions, and 

-4^"+' (m + 1)! (4n + m + l)!!!! (m-|)!!!! 



2(3 J (4n + m + 4)! (m+1)!!!! (4n + m - i)!!!! 
The notation a!!!! stands for: 



(a + 4)!!!! = (a + 4) .a!!!! . 
The convergence of the series, tested by the ratio test. 



lim 

n— »oo 





lim 






n- 


->oo 




Aji+i 






lim 


An 




2/9 







2/3 



lim 



0, 



' (An + m + T) 
{An + m + 2) (4n + m + 1) 



(4n + m + 5) (4n + m + 3) (4n + m - 2) 



0, 



^""^ (4n + m + 8) (4n + m + 6) (4n + m+ I) ' 
shows that both are convergent with an infinite radius of convergence. 

Therefore, in the first order approximation for /3, we have 



1 



1 



A (r) = -1^ + Ai + Asr + Asr^ + ^ (- 



(31) 



r 



6/3 V 4 10 60 175 
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An analysis of the solution (13T!) reveals the expected weak field behaviour at short scales 
{f), and deviat ion from this for a mesoscale, since we are dealing only with the linear 
approximation. Correspondingly, we find uq = 2GM/ where M is the mass of the central 
body. 

The remaining integration constants set scale distances where modifications of the Newto- 
nian behaviour appear. For instance, consider the Einstein-Hilbert theory with cosmological 
constant. The static spherically symmetric solution is 

, , ^ , , ^ 2GM 1 A 2 

v (r) = -A (r) = 1 — r 

r 6 

where the cosmological constant sets an scale distance given by the de Sitter pseudo-radius. 

Analogously, in our case, the vi constant sets a constant potential, which can be a mean 
nonlocal value of the effective Lagrangian proposed, V2 sets an scale distance where a constant 
mean force appears, and z/3 represents a gradient of force, in the same way as the cosmological 
constant in the example above. Similar reasoning can be developed for the other constants 
in the model. 

The contribution of each constant to the net force could be fixed by requiring that it 
fits the observational data for the tests of the gravitation. This task deserves a careful 
investigation of its own and is presently under investigation by the authors by means of the 
study of galaxy rotation curves, geodesic motion, perihelion shift, gravitational lenses and 
redshift. 



VII. CONCLUSION 



We have applied the second order gauge theory [10(] to the local gauge theory for the 
homogeneous Poincare group. It was found that the geometrical counterparts of the usual 
field strength F and the second order field strength G = DF are the Riemann tensor R 
and its (space-time) covariant derivative 5R. It followed the analysis of the second order 
invariants composed with the geometrical entities. 

We demonstrate - employing the symmetry properties of the curvature tensor - that the 
only independent Lagrangian densities for the gravitational field in a Riemannian manifold 
of arbitrary dimension are the seven ones listed in table (128|) . Linear combinations of terms 
proportional to powers of i?, as the familiar quadratic term in the curvature, are of first 
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order in the gauge potential cj, therefore, in the context of the second order gauge theory, 
the contributions of second order in the Lagrangian density, which are those including second 
derivatives of the gauge potential, are of type 5R. 

We derived equations of motion using a particularly simple choice for the second order 
gauge Lagrangian inspired in the Podolsky's proposal for a Generalized Electrodynamics. 
We found the static isotropic solution of these equations in the linear approximation, showing 
that at short distances the gravitational field behaves exactly as the Newton's law, but at 
meso-large distance scales the higher order contribution dominates, exhibiting a modified 
potential. 

In the future, we will study other solutions of these field equations, searching for massive 
modes which do not violate the local gauge symmetry. Our guide in these calculations 
shall be the treatment given in 1^ to the U (1) case, where an effective mass for the 
photon was derived. To do this, one naturally must concern about the determination of the 
conserved current associated with the local Lorentz symmetry and the relationship to the 
global diffeomorphic invariance of the theory. 

Another perspective is to apply the second order equations of motion ( l30l) to a Friedmann- 
Robert son- Walker metric. The goal is to seek for a accelerated regimes of the cosmological 
model arising from the higher order terms. This proposal is now under investigation. 
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APPENDIX A: APPENDIX: COUNTING SECOND ORDER INVARIANTS 
1. Counting GG Invariants 

First, let us analyse how many are the possible contractions of kind GG. This is done by 
means of tables as in the section IIV A[ The first one is constructed fixing, for instance, the 
last index: 



Fix. e 


Fix. a 


Fix. b 


Fix. c 


Fix. d 




(abed) 


(bacd) 


(cabd) 


(dabc) 


cyclic 


(acdb) 


{bcda) 


{cbda) 


{dbca) 




(adbc) 


(bdac) 


(cdab) 


(dcab) 




(abdc) 


(bade) 


{cadb) 


(dacb) 


non-cycl. 


(acbd) 


(bead) 


[chad) 


(dbac) 




(adcb) 


{bdca) 


{cdba) 


{dcba) 



Analogous tables result when we fix the indices d,c,b and a. For each table, non-cyclic 
permutations are equivalent to cyclic ones, giving: 



Fix. e 


Fix. a 


Fix. b 


Fix. c 


Fix. d 




[abed) 








cyclic 


(acdb) 


(bcda) 








(adbc) 


(bdac) 


{cdab) 





an similarly for the other four tables. 

Using the cyclic permutation symmetry, one can identify elements of different tables, 
reducing the number of invariants. By the skew-symmetry in the first G and renaming 
dummy indices, it follows: 



Qi = [abcde] [abode] 


Qq = [abcde] [cdbea] 


Q2 = [abcde] [beacd] 


= [abcde] [adbec] 


= [abcde] [adceb] 


Qs = [abcde] [acbde] 


= [abcde] [aecbd] 


Qq = [abcde] [acdeb] 


= [abcde] [debac] 


Qio = [abcde] [abdce] 



One can further apply the cyclic permutation symmetry to the first G in these remaining 
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invariants and reduce even more the number of independent quantities. Beginning with Qiq: 

Gio = — {[abdec] + [abecd]) [abdce] = Gi — Gio =^ ^Qiq = Qi . 

On the other hand, for Qg: 

Gg = — {[abdec] + [abecd]) [acdeb] = 2Q4 + Gg =^ = . 

Proceeding in the same way, one finds the foUowing identities: 

2^10 = Gi', 2^6 = G5 ', 

G2 = G3 = G4 = G7 = Gs = G9 = ^ ■ 

which give two independent invariants, 

if = [abcde] [abcde] , I2 = [abcde] [debac] . 



2. Counting (TrG) Invariants 

Starting with the three independent traces hsted in (!23|) . and considering skew- 
symmetries, the possible quadratic combinations are: 



TrGi = [-a ■ be] [-a ■ be] 


TrGu) = [ab ■ c] [ba ■ c] 


TrG2 = [■CL ■ be] [-b ■ ae] 


TrGii = [ab ■ c] [be ■ a] 


TrGi = [-ab ■ e\ [-a ■ be\ 


TrGi2 = [ab ■ e] [ea ■ b] 


TrGi = [-aft ■ c] [-b ■ ae] 


TrGis = [(^b ■ e] [eb ■ a] 


TrG^ = [-ab ■ c] [-c ■ ab] 


TrGu = [-ofe • c] [ab ■ -e] 


TrG 6 = [ab ■ -e] [-a ■ be] 


TrGi5 = [(^b ■ e] [ae ■ -b] 


TrGj = [ab ■ ■e\[-e ■ ab] 


TrGm = [o>b ■ c] [be ■ -a] 


TrG 8 = [ab ■ e] [ab ■ e] 


TrGn = [ab ■ -c] [ab ■ -c] 


TrGd = [ab ■ e] [ae ■ b] 


TrGis = [(lb ■ -e] [ae ■ -b] 



The last two invariants can not be converted into any other using the symmetries at our 
disposal. Each one of the preceding TrG must be analyzed, case by case, in a search for 
eventual interdependence. 

Take, for example, the 16**^ term, and rewrite it as bellow: 

TrGiQ = — [aeb-] [be ■ -a] — [a ■ eb] [be ■ -a] ^ 2TrGi6 = TrGr ■ 
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Repeat the reasoning for, say, the 15**^ invariant: 

TrQi5 = — [-acb-] [ac ■ -b] — [-a ■ cb] [ac ■ -b] = TrQu — TtQq . 

As soon as we perform this same check for all the above invariants, only eight of them 
are kept: 



Tr^3 = [ab ■ c] [a ■ be] 


TrQii — [-ab ■ c] [-be ■ a] 


TrQ^ = [-ab ■ c] [-c ■ ab] 


Tr^i4 = [-ab • c] [-ab ■ c] 


TrQa = [ab ■ ■c][-a ■ be] 


TrQn = [ab ■ -c] [ab ■ -c] 


TrQiQ — [-ab ■ c] [-ba ■ c] 


TrQis — [o^ • -c] [ac ■ -b] 



(Al) 



3. Counting (TrTrG) Invariants 



Prom T^ll = [-a ■ be] one can take a trace again: 



TW^[-o.oc] . 

Prom T^ll = [-ab ■ c] one finds Tc = [• o o • c] which can be reduced to T^^^ using the G 
skewsymmetry in the first two indexes and changing dummy indexes. Another possible trace 



is constructed from T 



(2). 



abc ■ 



.(2) 



[•o6-o] . 



But it also is not independent of tJ^^ : 



.(2) 



Let us set T^'^^ as the independent double trace. 



[• o 6 • o] = — [• o • o 6] — [• o o6- 



-T, 



(1) 



[0.0.6] 



-2T, 



(1) 



There are an internal double trace of T^H = [ab ■ -c] which is independent of Tc 



.(2). 



^ [ob . .0] 



The other double trace of Tjfl is. 



(3) 



n=[bo..o]^-T, 
Then, we have the following set of independent double traces: 

TrTrQi = [. o 6 . o] [. o 6 . o] 
TrTrQ2 = [°b ■ -o] [06 . .0] 
TrTrGs = [. o 6 . o] [06 . .0] . 



(A2) 



(A3) 



(A4) 
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4. Reducing the invariants using Bianchi identities 

Consider the reduction of the number of quadratic invariants in G by means of Bianchi 
identities. Using the geometric form, the first two invariants are: 

Applying the second Bianchi identity to we have: 

= -bfiR.p^^ [b'^R^'^P'' + b^R^'^f"') = If - If ^ 
^ 2lf = If , 

therefore it is sufficient to consider only If . 

Let us analyse now the trace invariants in G, (12^ : 



Trg^ = bfsR^^^ '^b,R"P>^^^ 


TrQu = bfsR^^b^R^^ 


Trg, = bpR^^^ Pb^RP^'^^ 


TrGu = bpR^J^R''^ 


TrGe = bPR^p^^b^R'^^^- 


TrGn = bPR^p^J^R"^^- 


Trg.o = bpR^^b^R^f" 


TrG^s = bPR^p^J^R^^^-- 



Comparing TrQiQ with TrQu one sees that both are the same invariant, due to the symmetry 
of Ricci tensor. 

Using the second Bianchi identity, it follows: 

Trg, = bpR^^g^p [b^RP^^^ + b'^R^P^'^) = TrQ,, - TrG.o 

and in the same way. 

This shows that only TrQi4 and TrQiQ can be hold independent. 
We apply the same technique to the double traced invariants f l25l) : 

TrTrGi = bpRbf^R 
TrTrG^ = bPRp^b^R^^ 
TrTrQs = -bpRb^R^^ . 
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The second Bianchi identity shows us that we have only one invariant in such case: 
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